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Section  1 


INTRODUCTION  AND  SUMMARY 


This  is  the  final  report  on  Contract  Number 
F19628-79-C-0125  titled  "Electron  Beam  Propagation  in  the 
Ionosphere" .  During  this  contract,  we  have: 

(1)  carried  out  a  fairly  thorough  literature 
survey  on  rocket  electron  (e)  beam  experi¬ 
ments  and  theoretical  investigations  dealing 
with  beam-plasma  interactions; 

(2)  developed  a  2-D  Monte  Carlo  electron  transport 
code; 

(3)  applied  the  code  to  beam  energy  deposition 
and  spreading;  and 

(4)  made  recommendations  with  regard  to  future 
experiments. 

The  motivation  behind  these  efforts  is  an  Air  Force  Geophysics 
Laboratory  program  designed  to  fire  kilovolt  ion  and  electron 
beams  from  rockets  and  possibly  the  space  shuttle. 

The  highlight  of  the  work  performed  was  the  develop¬ 
ment  and  application  of  code  MCBE  which  is  the  code  referred 
to  above.  This  code  solves  the  equation  of  motion  for  some 
large  number  of  electrons  in  the  presence  of  an  exponentially 
varying  neutral  atmosphere  and  geometric  field.  Along  the 
electron's  path  of  motion,  energy  is  allowed  to  be  deposited 
continuously.  Frequent  changes  in  the  electrons'  guiding 
center  take  place  through  scattering  which  is  modeled  by 
use  of  the  Goudsmit-Saunderson  multiple  scattering  formula. 
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The  Monte  Carlo  nature  of  MCBE  comes  from  scattering 
treatment  and  in  some  cases  from  sampling  distribution  functions 
defining  the  initial  beam  conditions.  The  quantities  calculated 
are: 

(1)  the  r,z  distributions  of  energy  and  charge 
deposition; 

(2)  the  1-D  distributions  of  energy  deposition 

in  either  r  or  z  obtained  by  integrating  over 
the  other  spatial  dimension;  and 

(3)  detailed  radial  and  time  information  on  those 
electrons  which  have  backscattered  through 
the  incident  plane. 

The  current  description  of  electron  interactions  in  MCBE 
does  hot  include  wave-particle  effects. 

A  noteworthy  result  obtained  from  various  runs 

of  MCBE  is  the  insensitivity  of  the  spreading  to  the  starting 

energy  of  the  beam  electrons.  This  contrasts  with  the  results 

24 

of  Berger,  et  al .  (see  the  appendix  for  further  details) 
and  leads  to  differences  on  the  order  of  three  in  mean  radius 
of  energy  deposition  for  a  starting  energy  of  20  keV.  Because 
of  this  disagreement,  an  analytic  model  of  guiding  center 
diffusion  was  developed  to  test  the  MCBE  results.  The  model 
confirms  both  the  insensitivity  of  spreading  to  the  starting 
beam  electron  energy  and  the  magnitude  of  the  spreading. 

This  work  has  been  published  with  a  copy  of  the  reprint  appearing 
in  the  appendix. 

An  application  of  MCBE  and  the  above  analytic 
model  has  been  made  to  interpret  the  observed  optical 
streak  widths  produced  by  ECHO-4  beams  emitted  downward  into 
the  E-region.  The  calculated  widths  give  reasonable  agreement 
with  the  observations  for  some  streaks  and  significantly 


underestimate  the  observations  for  others.  Generally,  beam 
spread  increased  for  those  beams  fired  near-rocket  apogee. 

This  suggests  that  the  differences  can  be  attributed  to  plasma 
turbulence  effects  occurring  in  regions  of  lower  collision 
frequency.  The  comparisons  quantify  the  amount  of  the  beam 
spreading  due  to  the  turbulence. 

Each  of  the  chapters  to  follow  was  originally  written 
as  a  self-contained  document.  Section  2  comes  from  an  interim 
report  which  documents  our  literature  search.  Section  3 
comes  from  combined  Quarterly  Reports  2  and  3  and  presents 
MCBE  results  and  relates  them  to  ECHO-4  data.  A  discussion 
is  included  of  the  rate  of  plasma  turbulence  in  the  spreading 
of  the  emitted  beams.  Section  4  contains  a  memo  submitted 
during  the  contract  period  providing  suggestions  for  future 
rocket  fixed  beams  experiments.  Finally,  the  appendix  contains 
a  copy  of  a  JGR  paper  based  on  part  of  the  contract  work. 

The  paper  describes  the  analytic  model  for  beam  spreading 
referred  to  above. 
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Section  2 

A  REVIEW  OF  ROCKET  e-BEAM  EXPERIMENTS 
AND  ASSOCIATED  THEORETICAL  INVESTIGATIONS 

2.1  INTRODUCTION 

Since  the  first  rocket  experiment  with  injection 
of  e-beams  into  the  ionosphere  by  Hess  et  al.^\  nearly 
two  dozen  such  experiments  have  been  conducted,  providing 
an  impressive  body  of  information  on  e-beam  propagation 
in  the  ionosphere  and  magnetosphere.  This  information  has 
come  from  combinations  of  on-board,  separated  nose  cose, 
and  ground  based  measurements  of  particles,  optical  emissions, 
and  rf  emissions.  Much  has  been  learned  about  the  interaction 
of  e-beams  with  the  neutral  particles  and  plasma  within 
the  propagating  region,  but  such  experiments  have  also  left 
a  number  of  questions  directly  related  to  beam-plasma  inter¬ 
actions  unanswered. 

We  will  attempt  in  this  report  sponsored  by  the 
Air  Force  Systems  Command  under  Contract  Number 
F19628-79-C-0125 ,  to  provide  a  fairly  comprehensive  review 
of  the  above  subject,  noting,  by  experiment,  parameters 
of  interest  and  those  observed  phenomena  which  are  not  cur¬ 
rently  well  understood. 

2.2  ROCKET  e-BEAM  EXPERIMENTS 

2.2.1  Overview 


To  date  (late  1979)  we  are  aware  of  fifteen  rocket  e-beam 

experiments  that  have  been  conducted  in  the  ionosphere. 

These  are  two  early  experiment  reported  by  Hess  et  al.^\ 

(  2 )  (  3 ) 

Davis  et  al.  ,  and  Davis  ,  four  experiments  labeled 
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ECHO  I-ECHO  IV,  the  Russian  experiments  ZARNITSA  1-2,  the 
Franco-Russian  experiment  ARAKS,  five  experiments  in  the 

Q  } 

PRECEDE-EXCEDE  program,  and  a  Norwegian  experiment  .  In 

(5) 

addition,  Cohen  et  al .  have  conducted  an  important  experi¬ 
ment  where  a  rocket  fired  ion  beam  (in  conjunction  with 
an  electron  beam)  was  used  to  study  the  effects  of  beam 
emission  on  spacecraft  charging.  Since  characterization 
of  beam  propagation  was  not  a  primary  goal  of  the  experiment, 
it  has  not  been  considered  in  this  review.  The  value  of 
beam  propagation  experiments  lies  in  their  potential  to 
obtain  information  on  such  subjects  as: 


•  the  geomagnetic  field: 

conjugate  point  locations, 
field  line  lengths, 
magnetospheric  field  configuration 
(e.g.,  transition  from  closed  to 
open  field  lines); 

•  drift  rates: 

curvature, 
jjjrad^ent , 

•  particle-particle  interactions: 

transport  properties  such  as  scattering 

and  penetration  depths, 

neutral  densities, 

cross  sections, 

rate  coefficients; 

•  wave-particle  interactions. 

The  fundamental  problem  of  e-beam  propagation 
comes  under  the  last  two  headings  and  addresses  how  particle- 
particle  and  wave-particle  interactions  modify  the  beam's 


( 4 ) 

a)  Recently,  Winckler  has  written  a  report  which  expands 
this  list  by  including  the  Russian  "Feyerwerk"  experiment, 

5  experiments  in  the  Japan  "K"  rocket  series,  and  2  E||B 
experiments  with  rockets  fired  from  the  Churchill  Research 
Range.  Information  on  the  experiments  is  not  now  available 
in  the  open  literature,  and  they  have  not  been  considered 
in  this  report. 
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lateral  extent  and  its  distributions  in  energy  and  pitch 
angle.  Prior  to  the  first  e-beam  experiment  in  the  iono¬ 
sphere  (Hess  et  al.^^),  it  was  an  open  question  as  to  whether 
a  beam  would  even  propagate  an  appreciable  distance  before 
some  plasma  instability  led  to  its  destruction.  Experiment 
has  shown  that,  indeed,  beams  do  propagate  and  seem  to  retain 
much  of  their  character  over  great  distances  such  as  traversed 
in  an  "echo"  from  the  conjugate  point.  Experiment  has  also 
shown  that  mechanisms  are  operating  on  the  beam,  presumably 
of  the  collective  or  wave-particle  type  which  lead  to  effects 
poorly  understood  at  this  time.  Examples  are  anomalous 
backscatter,  high  frequency  rf  generation,  and  anomalous 
diffuseness  of  artificial  auroral  streaks. 

In  the  sub-sections  to  follow,  we  will  be  addressin 
the  above-mentioned  experiments  by  way  of  presenting  launch 
information,  beam  parameter  values,  diagnostics,  and  unexplained 
effects.  We  do  wish  to  caution  the  reader  that  our  information 
comes  from  published  literature  and  is  not  necessarily  complete. 
An  example  of  this  might  be  pulse  durations  where  in  fact, 
the  values  given  may  be  part  of  a  larger  set  not  discussed 
in  its  entirety  in  the  literature.  We  believe,  however, 
that  sufficient  information  will  follow  to  provide  the  reader 
with  a  good  picture  of  rocket  e-beam  experiments  and  what 
is  and  is  not  understood  about  them  to  date. 

2.2.2  Parameters  and  Instrumentation 


An  attempt  has  been  made  by  us  to  characterize 
the  various  experiments  in  the  form  of  tables  giving  a  variety 
of  paramaters  and  instruments  which  recorded  data  of  a  general 
interest  to  this  program.  The  reader  will  find,  however, 
a  number  of  blank  spaces  in  these  tables,  especially  for 
instrumentation.  This  simply  means  that  we  have  not  found 
the  given  information  in  the  literature  base  used  for  this 
review.  An  advantage  of  the  format  presentation  chosen. 
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in  spite  of  this,  is  that  newly  acquired  information  may 
be  easily  added  which  we  intend  to  do  as  this  program  continues 

Table  1  provides  a  list  of  the  experiments  along 
with  basic  information  such  as  launch  date,  launch  site, 
and  rocket  apogee.  We  see  that  the  experiments  cover  a 
time  span  of  one  decade  and  latitude  range  from  mid  to  high 
latitude.  This  latter  fact  is  worth  noting  since  environmental 
conditions,  in  terms  of  beam-plasma  interactions,  may  vary 
significantly  from  the  mid-latitude  ionosphere  to  the  auroral 
ionosphere  where  natural  streams  of  energetic  electrons 
are  frequently  present. 

The  parameters  we  have  chosen  to  characterize 
the  beam  are  energy,  current,  power,  pulse  duration,  and 
injection  angle.  An  example  of  how  some  of  these  parameters 
have  been  programmed  into  an  experiment  (ECHO  IV)  is  provided 
in  Figure  1.  We  observe  a  range  of  pulse  widths  and  injection 
pitch  angles  which  is  not  atypical  for  rocket  e-beam  experi¬ 
ments.  More  information  is  given  in  Table  2  where  we  have 
attempted  to  specify  the  above  parameters  by  experiment. 

We  wish  to  stress  that  the  parameter  values  shown  may  be 
parts  of  larger  sets.  We  believe,  however,  that  most  entries 
do  cover  the  experimental  range.  The  convention  chosen 
on  pitch  angle  is  that  0°  corresponds  to  the  downward  direction 

The  basic  types  of  recording  devices  used  in  the 
experiments  were: 

•  particle  detectors, 

•  photometers, 

•  photo-spectrometers, 

•  rf  receivers, 

•  radar  systems,  and 

•  TV  systems. 
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TABLE  1 


ROCKET  e-BEAM  EXPERIMENTS 


LAUNCH  AND 


TRAJECTORY  PARAMETERS 


HESS  ET  AL 
DAVIS 
ECHO  I 
ECHO  II 
ECHO  III 
ECHO  IV 
ZAHNITSA-I 
ZARNITSA-II 
ARAKS 

(2  experiments) 


BK3DE-II  TEST 


PRECHK-II 

EXCEDE-SPECITtA 

NORWEGIAN 


LAUNCH 

DATE 

LAUNCH 

SITE 

L 

PARAMETER 

Jan.  '69 

Wallops  Is.,  Va. 

Kauai,  Hawaii 

2.6 

Aug. 

•70 

Wallops  Is. ,  Va. 

2.6 

Sept.  '72 

Ft.  Churchill,  Manatoba 

•v8. 8 

Apr. 

'74 

Poker  Flat,  Alaska 

6 

Jan. 

•76 

Poker  Flat,  Alaska 

6 

Jan. 

'75 

Kerguelen  Island 

Feb. 

•75 

(49*  S) 

Oct. 

'74 

White  Sands,  N.  M. 

•vB 

Apr. 

•75 

Poker  Flat,  Alaska 

6 

Feb. 

•76 

Poker  Flat,  Alaska 

6 

Dec. 

'77 

White  Sands,  N.  M. 

^ 2 

Oct. 

•79 

Poker  Flat,  Alaska 

6 

FIGURE  1.  ECHO  IV  Gun  Sequence 
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TABLE  2.  ROCKET  e-BEAM  EXPERIMENTS  —  BEAM  PARAMETERS 


EXPERIMENT 

MAX.  POWER 
(kW) 

PULSE 

DURATION 

(ms) 

BEAM 

ENERGY 

(keV) 

BEAM 

CURRENT 

(tnA) 

INJECTION 

ANGLES 

(DEGREES) 

HESS  ETT  AL 

5 

<1000 

<10 

<500 

DAVIS 

ECHO  I 

3 

16 

35-43 

70 

68-112 

EDO  II 

3.8 

64 

4-47 

80 

ME- 180 

ECHO  III 

2,8,32 

32-37 

ECHO  IV 

2.5 

1000 

16-38 

70 

70,110,30-150 

(sweep-node) 

ZAHNITSA-I 

ZARNITSA-II 

4.5 

400 

7,9 

300-500 

30-90 

ARAKS 

13.5 

20,1280, 

2560 

15,27 

<500 

0,70,140 

JW3CEXJE 

2 

2.5 

800 

140-180 

BtCEDE-II  TEST 

30 

M.000 

EXCEDE-SWIR 

3 

3 

1000 

PRECEEE-II 

20 

4300 

3 

7000 

EXCEDE-SPECTRA 

NORWEGIAN  EXPERIMENT 

In  a  given  experiment,  some  of  these  would  be  located  on 
the  ground  around  the  launch  site,  on  the  payload,  and  possibly 
on  the  ground  below  the  conjugate  point.  Table  3  provides 
a  list  of  instruments  by  experiment  and  places  them  into 
either  an  on-board  or  ground-based  category. 

Several  types  of  particle  detectors  are  noted 
in  the  table.  Most  of  these  were  used  to  measure  the  return 
of  beam  particles  due  to  mirroring  and  scattering  (by  particle- 
particle  and  wave-particle  interactions).  Most  of  the  effort 
which  went  into  the  ECHO  experiments  was  directed  to  this 
type  of  measurement.  The  retarding  potential  analyzer  (RPA), 
however,  was  used  to  measure  the  return  current  (low  energy) 
needed  to  keep  the  electron-guns  functioning.  On  all  experi¬ 
ments,  attention  was  given  to  vehicle  potential  and  its 
variation  with  time. 

RF  recievers  were  used  to  record  signals  generated 
by  the  beam  interacting  with  the  plasma.  In  some  cases, 
these  receivers  were  well  positioned  by  being  located  in 
ejected  nose  cones  (e.g.,  on  ECHO  I,  ECHO  II,  and  ARAKS). 

For  those  cases,  rf  signals  able  to  propagate  in  the  plasma 
environment  of  the  ionosphere  could  be  detected.  Observations 
of  ground  based  rf  signals  have  been  more  limited,  but  radio 
signal  signatures  from  beam  injection  experiments  have  been 
obtained. 


The  remaining  instrumentation  appearing  in  Table  3 
needs  litte  explanation.  Radar  systems  were  used  to  track 
the  rocket  (not  relevant  to  this  discussion)  and  measure 
ionization  produced  by  beam  energy  degradation  in  the  lower 
ionosphere  (e.g.,  in  the  ARAKS  experiment).  Photometers, 
usually  recording  Ng  3914A  and  01  5577A  emission,  provided 
a  measure  of  how  much  energy  was  being  deposited  when  the 
beam  produced  an  artificial  auroral  streak.  Photo-spectro¬ 
meters  give  the  strengths  of  many  optical  features  which 
may  be  used  to  better  understand  chemical  processes.  This, 
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TABLE  3  ROCKET  e-BEAM  EXPERIMENTS— INSTRUMENTATION 


EXPERIMENT 

ON-BOARD 

GRCKJND-RASED 

— 

Retarding  Potential  Analyser  (RPA)  (0-2  keV) 

TV,  Radar,  RF,  Photometers 

DAVIS 

TV 

ECHO- 1 

RF  (in  ejected  nose  cone) 

Particle  (scintillation,  solid  state) 

ECHO- 1 1 

RF  (in  ejected  nose  cone), 

Particle  (electrostatic  (.9-45  keV),  solid  state) 

Proportional  (for  Bremsstrahlung) 

ECHO-III 

Particle  (electrostatic,  scintillation) 

Photometer  (3914A) 

ECHO  IV 

Particle  (unspecif  ied\Haotoroeters  (3914A  &  5577A) 

TV 

ZARNITSA-I 

Particle  (electrostatic  l-l4  keV) , Magnetometer, 
Radar  Transponder 

TV,  Radar, RF( 300Hz-15KHz) 

ZARNITSA-II 

TV,  Radar,  RF  (22.5-80MHz) 

ARAKS 

Particle,  RF  (in  ejected  nose  cone) 

Radar  (conjugate  point) 

RF  (10-75  MHz) 

PRECEDE 

RPA 

TV,  Photometers  (3914A  &  5577A) 

Spectograph (4200-8500A) 

EXCEDE-II 

TEST 

RPA,  Photometers  (3914A  &  5577A), 

Radiometers  (2.7  &  4.3  \m) 

TV,  Photometers  (3914A&  5577A) 

EXCEDE-SWIR 

22  Photometers,  Radiometer,  Ion  Mass  Spectrometer 

TV,  Photometers  (3914A&  5577A) 

PHECEDE-II 

IR  Interferometer 

TV,  Ihotemeters  (3914A  &  5577A) 

UV  &  Visible  Spectrographs 

EXCEDE- 

SPECTRA 

in  fact ,  has  been  the  emphasis  in  the  PRECEDE-EXCEDE  program 
with  particular  attention  directed  to  chemistry  affecting 
IR  emission.  Finally,  TV  systems  have  been  employed  to 
search  for  auroral  streaks  produced  by  either  downward  injec 
tion,  upward  injection  from  low  in  the  ionosphere,  or  upward 
injection  at  higher  altitudes  leading  to  an  echo.  Streaks 
have  been  recorded  on  many  of  the  experiments  starting  with 
the  first  one  reported  by  Hess  et  al.^^  In  none  of  the 
experiments,  however,  have  streaks  been  recorded  for  beam 
energy  deposition  following  an  echo. 

2.2.3  Observations 


Based  on  the  literature  base  we  have  developed 
to  this  point ,  the  experiments  of  most  interest  to  this 
program  are  ECHO  I  -  IV ,  ZARNITSA-2  and  ARAKS.  Taken  as 
a  whole,  these  experiments  (as  well  as  the  first  by  Hess) 
have  significantly  advanced  our  understanding  of  e-beam 
propagation  in  the  ionosphere  and  magnetosphere.  They  have 
also  obtained  several  results  which  cannot  be  explained 
by  classical  theory.  We  wish  to  point  out  such  results 
in  this  section  and  to  do  so  we  have  constructed  Table  4. 
Here,  we  have  divided  the  observations  into  the  categories 
of  particle,  rf,  and  optical  for  six  experiments  noted  above 
Blank  spaces  mean  either  the  absence  of  that  particular 
measurement  or  the  absence  of  information  in  the  literature. 
In  some  places  we  have  noted  the  measurement  even  if  it 
produced  no  surprises.  Our  emphasis  here,  however,  is  on 
measurements  which  imply  strong  beam-plasma  interactions 
and  thus  in  turn  are  not  well  understood.  Prior  to  their 
discussion,  however,  we  note  that  particle  measurements 
such  as  those  on  ECHO  IV  show  that  the  beam  can  remain 
reasonably  well  confined  over  large  propagation  distances 
such  as  traversed  in  an  echo  from  the  conjugate  region. 

These  imply  that  beam-plasma  interactions  are  not  particu¬ 
larly  destructive  to  the  beam.  There  is  reason  to  believe. 


Table  4.  ROCKET  e-BEAM  EXPERIMENTS  -  OBSERVATION  HIGHLIGHTS 


however,  that  the  strength  of  the  interaction  can  vary  sig¬ 
nificantly  depending  on  the  specific  properties  of  the  propa¬ 
gation  medium  and  particularly  on  the  position  of  beam  injection. 

A  feature  requiring  the  presence  of  a  rather  strong 
beam-plasma  interaction  is  backscatter  from  just  above  the 
rocket  as  observed  on  ECHO  I  and  II  and  on  ARAKS.  The  measured 
signals  were  orders  of  magnitude  above  that  predicted  for 
particle-particle  scattering.  Another  such  feature  is  the 
halo  of  strong  enhanced  ionization  (not  accounted  for  by 
injected  plasma  sources)  in  the  rocket  environment  as  observed 
on  the  ZARNITSA  and  ARAKS  experiments.  The  halos  were  detected 
by  high  frequency  rf  signals  propagating  from  the  vicinity 
of  the  rockets  and  by  ground  based  radar  and  TV  systems. 

In  the  ZARNITSA  2  experiment,  it  was  suggested  that  several 
percent  of  the  initial  beam  energy  became  converted  to  plasma 
potential  and  kinetic  energy  in  the  local  environment  of 
the  rocket.  A  number  of  other  observations  can  be  cited 
supporting  significant  beam-plasma  interactions.  Diffuse 
auroral  streaks  were  observed  on  ECHO  IV  when  the  beam  origi¬ 
nated  near  apogee  (215  km).  Unexpected  delays  in  echos 
were  observed  on  ECHO  I  and  III.  In  the  ARAKS  experiment 
long  delays  in  the  onset  of  radar  return  signals  from  regions 
of  beam  deposition  at  the  conjugate  point  were  observed. 

In  some  instances  of  ECHO  I  injection  there  were  observed 
pairs  of  echos  with  the  first  arriving  at  about  the  predicted 
time  while  the  second  arrived  much  later  but  too  soon  to 
be  from  a  second  conjugate  bounce. 

In  addition  to  observations  that  show  beam  dissi¬ 
pation,  rf  signals  have  been  observed  to  result  from  beam- 
plasma  interactions  which  largely  maintain  the  integrity 
of  the  beam.  These  signals,  although  not  now  well  understood, 
could  prove  to  be  valuable  diagnostics  on  the  characteristics 
of  the  propagating  beam.  The  frequencies  of  these  signals 
are  in  the  MHz  range  corresponding  to  ambient  plasma  frequencies 
as  well  as  electron  cyclotron  harmonics  and  in  the  sub- 
MHZ  range  corresponding  to  whistler  modes. 
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2.3 


COLLECTIVE  EFFECTS 


2.3.1  Overview 

By  collective  effects  we  mean  the  totality  of 
fluctuations  in  a  plasma  with  which  a  rocket  fired  beam  may 
interact  to  result  in  processes  different  than  or  substantially 
enhanced  from  usual  particle-particle  interactions.  Since 
the  ambient  properties  of  the  ionosphere  and  magnetosphere 
can  be  determined  reasonably  well,  it  has  been  possible 
to  anticipate  the  results  of  particle-particle  interactions 
as  well  as  the  classical  drift  character  of  the  individual 
beam-particles  within  inhomogeneous  magnetic  fields.  It 
is  the  complex  array  of  unanticipated  results  which  have 
occurred  in  rocket  fired  beam  experiments  that  has  given 
rise  to  an  intense  interest  in  collective  effects. 

In  the  study  of  these  effects  it  is  useful  to 
consider  three  different  regions  of  space.  The  first  is 
the  region  close  to  the  beam  injector  where  potential  build¬ 
up  between  various  rocket  components  is  possible.  Such 
configurations  can  lead  to  important  return  current  effects, 
which  often  can  be  enhanced  by  the  presence  of  rocket  gases 
remaining  in  the  vicinity  of  the  gun.  Not  only  is  it  difficult 
to  characterize  the  environment  of  the  gun  but  it  can  be 
expected  that  disruptive  beam-particle  interactions  will 
occur  near  the  injection  position.  The  major  phenomenon 
observed  near  the  rocket  is  an  enhancement  of  the  ambient 
electron  density  in  an  effect  known  as  the  beam  plasma  discharge 
(BPD).  The  second  region  of  interest  is  the  ionosphere 
at  some  distance  from  the  rocket.  Here  the  ambient  environment 
is  often  times  sufficiently  free  of  fluctuations  that  we 
can  confidently  associate  collective  effects  with  fluctuations 
induced  by  the  beam-plasma  interaction.  The  effects  of 
beam  interaction  within  the  ionosphere  have  been  monitored 
through  nose  cone  antennas  which  receive  rf  signals.  Among 
suspected  collective  effects  have  been  the  appearance  of 
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diffuse  artificial  aurora  streaks  as  well  as  enhanced  local 
backscatter  with  either  upward  or  downward  injection.  The 
third  region  of  interest  is  the  magnetosphere  through  which 
the  beam  propagates  in  experiments  involving  echos  or  artificial 
aurora  created  at  the  conjugate  point.  Collective  fluctuations 
in  this  region  may  be  caused  by  instability  induced  through 
natural  phenomena  as  well  as  through  beam-plasma  interaction. 

The  major  observation  of  collective  effects  in  this  region 
involve  echo  delays  as  well  as  large  degradations  in  the 
expected  flux  of  mirroring  particles. 

In  Section  2.3.2  we  examine  the  important  parameter 
regimes  in  space.  The  importance  in  differentiating  between 
absolute  and  convective  instabilities  in  finite  width  beam 
experiments  is  discussed  in  Section  2.3.3.  Classification 
of  the  rf  signals  received  by  antennas  with  respect  to  plasma 
modes  is  made  in  Section  2.3.4.  Finally,  in  Section  2.3.5 
we  discuss  possible  mechanisms  to  explain  the  experimental 
results  apparently  involving  collective  effects.  Unfortunately, 
no  concensus  seems  to  exist  as  to  the  cause  of  the  major 
collective  phenomena.  This  points  up  the  importance  not 
only  of  increasing  the  emphasis  on  theoretical  studies  but 
also  in  the  design  of  more  controlled  experiments  in  the 
future . 

2.3.2  Parameter  Regimes 

When  one  considers  the  nature  of  beam-plasma 
collective  effects  that  can  arise  in  experiments  involving 
beams  fired  from  rockets,  it  is  important  to  keep  in  mind 
the  wide  range  of  parameter  regimes  that  are  encountered. 

Beams  are  generally  ejected  from  rockets  between  the  E  and 
F  regions.  Table  5  indicates  the  neutral  and  electron  number 
densities  at  the  various  heights  when  beams  are  injected. 

At  low  injection  altitudes,  the  neutral  particle  density 
may  be  so  large  that  a  very  high  neutral -electron  collision 
frequency  may  be  thought  to  preclude  any  collective  effects. 
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TABLE  5.  TYPICAL  IONOSPHERIC  PARAMETERS  AS  FUNCTIONS  OF  HEIGHT  (L=2.56) 


This  can  be  misleading,  however.  Even  though  the  high  collision 
frequency  may  serve  to  prevent  the  growth  of  large  amplitude 
waves  (ignoring  for  the  moment  resistive  instability),  still, 
self-electric  field  effects  can  be  important,  since  the 
background  plasma  density  decreases  for  lowering  altitudes 
beneath  the  E-layer.  The  usual  test  for  the  importance 
of  the  self-electric  field  is  to  compare  the  beam  density 
with  the  ambient  plasma  density.  The  radius  of  the  beam 
in  the  vicinity  of  the  rockets  can  usually  be  assumed  to 
be  or  order  of  the  Larmor  radius: 
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where  V  is  the  speed  of  an  ipjected  beam  particle,  0  is 
the  pitch  angle,  and  w  is  the  electron  gyro-frequency. 

Since  for  the  non-relativistic  particle  beams  of  usual  interest 
the  speed  is  related  to  the  energy  by  : 

v  -VF~ 


where  M  is  the  electron  mass,  and  since  the  beam  number 
density  is  related  to  the  total  current  by: 

2 

N^eVcosO  •  ffP  =  I  , 


we  can  finally  express  Nb  in  terms  of  the  usual  parameters 
associated  with  the  experiment  as  well  as  with  the  ambient 
environment : 
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Equation  (2)  shows  that  for  the  practical  case  of  small 
pitch  angle  will  be  larger  than  the  ambient  plasma  density. 

In  this  case  the  ambient  plasma  cannot  charge  neutralize 

the  beam  so  that  the  electrostatic  repulsion  should  be  included. 
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Gendrin'  and  Alekkin  et  al .  have  examine  the  radial 
spread  by  order  of  magnitude  arguments  and  conclude  that 
beam  spreading  ceases  after  the  perpendicular  velocity  attains 
a  certain  (small)  value.  The  final  average  number  density 
in  the  beam  is  shown  by  Gendrinv  '  to  correspond  to  a  plasma 
frequency  about  equal  to  the  electron  cyclotron  frequency. 

One  of  the  key  parameters  determining  the  behavior 

of  collective  effects  is  the  ratio  of  the  plasma  frequency 

to  the  electron  cyclotron  frequency.  For  altitudes  above 

(  8  ) 

the  E-layer  considered  by  Jones  and  Kellogg  ,  this  ratio 
is  generally  greater  than  unity  and  approaches  a  maximum 
value  in  the  F-layer.  In  the  E-layer  itself  the  ratio  is 
much  larger  than  unity.  For  altitudes  beneath  the  E-layer, 
however,  the  ratio  decreases  sharply  as  the  plasma  density 
decreases.  The  low  altitude  region  where  a)pe/a)ce<l  was 
of  considerable  interest  in  the  low  trajectory  flight  of 
ZARNITSA-2.  As  we  will  show  below,  this  is  the  same  parameter 
regime  which  describes  the  recent  series  of  beam  plasma 
discharge  (BPD)  experiments  carried  out  at  Johnson  Space 
Flight  Center. 

In  Table  5  we  present  data  on  ambient  conditions 
along  a  typical  magnetic  field  line  extending  from  the  low 
ionosphere  up  through  the  magnetosphere.  This  table  will 
be  a  useful  reference  in  the  determination  of  the  importance 
of  collective  effects  for  experiments  of  interest. 

2.3.3  Absolute  and  Convective  Instabilities 

The  behavior  of  the  rocket  fired  electron  beam 
depends,  of  course,  on  the  environment  through  which  it 
passes.  In  an  Echo  experiment,  for  example,  where  the  beam 
propagates  through  the  magnetosphere,  the  beam  may  encounter 
regions  of  turbulence  established  by  naturally  occurring 
streams  of  particles.  Such  turbulent  regions  may  scatter 
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as  well  as  change  the  energy  of  the  beam  particles.  The 
effect  of  such  processes  depends  on  the  origin  of  the  turbulence. 
It  is  for  this  reason  that  the  study  of  beam  interaction 
with  the  magnetospheric  plasma  requires  a  knowledge  of  the 
multitude  of  natural  phenomena  occurring  in  that  region. 


For  the  case  that  the  environment  can  be  considered 
relatively  stable,  the  beam  can  create  its  own  level  of 
turbulence  with  which  it  can  interact  non-linearly .  Unlike 
the  broad  streams  of  plasma  associated  natural  phenomena, 
rocket  fired  beams  have  relatively  small  diameters,  and 
this  feature  of  radial  inhomogeneity  can  be  important  in 
determining  whether  the  beam  plasma  interaction  will  drive 
parts  of  the  wave  spectrum  unstable.  The  radial  extent 
of  the  beam  usually  can  be  taken  to  be  of  order  of  the  electron 
Larmor  radius,  and  therefore  a  mode  propagating  in  the  beam 
plasma  system  must  have  a  radial  wavenumber  at  least  as 
large  as  the  inverse  radius.  In  many  cases  this  constraint 
on  the  perpendicular  wavenumber  can  have  a  stabilizing  effect. 
The  finite  radius  also  means  that  it  is  possible  for  an 
unstable  wavepacket  to  propagate  across  the  system  before 
the  instability  has  a  chance  to  grow  substantially.  Thus, 
when  the  instability  is  able  to  convect  out  of  the  region 
of  instability,  the  level  of  turbulence  will  remain  very 
low.  On  the  other  hand,  the  wavepacket  will  be  able  to 
grow  appreciably  when: 


>>  1 


where  pj_  is  the  beam  radius,  y  is  the  growth  rate  of  the 
mode,  and  V  is  the  perpendicular  group  velocity.  F  is 
a  measure  of  the  gain  in  the  amplitude  of  the  disturbance 
as  it  propagates  across  the  beam.  In  the  limit  of  V  -K> 
a  wavepacket  does  not  convect  away  no  matter  how  small  the 
region  of  instability  or  the  growth  rate  is.  Such  an  insta¬ 
bility  has  been  termed  absolute  (as  opposed  to  convective 
for  non-zero  Vg)  by  Briggs^®\  who  did  pioneering  work  in 


this  field  and  pointed  out  a  graphical  prescription  for 
determining  when  Vg  would  vanish.  Briggs'  representative 
construction  occurs  as  a  consequence  of  the  inverse  Fourier 
transform  used  to  obtain  the  spatial  distribution  of  an 
unstable  wavepacket  in  the  perpendicular  (x)  direction. 

If  a  change  of  variables  is  made  so  that  the  integral  over 
Kx  goes  into  an  integral  over  u>,  then  the  Jacobian  factor 

3K  (o))/3 to  is  introduced  into  the  integral.  The  factor  is 

0) 

formed  through  the  requirement  that  the  dielectric  constant 
e  (K  ,K  ,cj)  remains  zero  as  w  and  K  are  varied: 

O  X  Z  X 

eo<Kx’Kz’w)  ^  1 nr  <“-*>  +  TT  (Kx  -  *x> 

+  higher  order  terms  =  0 
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where  u>  and  K  =  K  X  +  K  Z  are  the  frequency  and  wavenumber 

X  z 

of  a  growing  mode  (lmuj>0  for  a  mode  with  exp(-iwt)  time 

/\ 

dependence),  and  Z  is  the  unit  vector  in  the  magnetic  field 
direction.  For  this  case  we  have  the  following  single  valued 
relation  for  K  as  a  function  of  u>: 

X 

K  -  g  =  (u-u)/qr 
x  x  ugj_ 


and  the  group  velocity  V  =  R  for  )•  This  type  of  relation 

®  ©JL 

for  a  non-zero  group  velocity  is  characteristic  of  a  convec¬ 
tive  mode.  When  or  -*-0  so  that  V„  also  -K),  however,  a  dif- 

ug±  g_L 

ferent  relation  between  Kx  and  u>  results.  For  this  case: 
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+  higher  order  terms  =  0 
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so  that  now: 


so  that  there  now  exist  two  branches  of  Kx(w)  in  the  vicinity 
of  w  and  K  .  The  above  behavior  holds  for  w  and  K  very 

X  X 

close  to  the  unstable  mode  characterized  by  u)  and  Kx-  The 
method  of  Briggs  generalizes  the  above  approach  by  examining 
the  behavior  of  Kx  over  a  prescribed  region  of  to.  For  the 
case  that  3to/9Kx  =  0  within  the  interval  of  interest,  a 
saddle-point  will  indicate  the  presence  of  the  absolute 
instability . 

The  presence  of  convective  rather  than  absolute 

instabilities  in  finite  radius  beam-plasma  interaction  allows 

for  an  effective  stabilizing  mechanism.  Indeed,  Jones  and 
(  8  ) 

Kellogg  used  the  behavior  of  convective  unstable  whistler 
modes  to  argue  that  beams  in  the  upper  ionosphere /lower 
magnetosphere  could  retain  their  integrity  while  propagating. 
The  analysis  presented  above  also  can  be  used  to  demonstrate 
the  possible  existence  of  steady  state  configurations  where 
the  growth  is  spatial  rather  than  temporal.  This  situation 
exists  when  the  dispersion  relation  is  solved  for  K  in 
terms  of  real  u)  rather  than  the  other  way  around.  It  is 
clear  that  finite  spatial  growth  is  only  meaningful  for 
convective  modes,  where  a  unique  relation  for  2  in  terms 
of  w  exists.  Note  that  the  relation  between  the  temporal 
growth  rate  y  and  the  spatial  growth  is  just 
y  =  so  that  the  concept  of  spatial  growth  is  not 

valid  for  absolute  instability  when  =  0.  The  idea  of 
steady  state  turbulent  levels  described  by  mode  amplitudes 
obtained  from  spatial  growth  considerations  in  the  inhomoge¬ 
neous  magnetosphere  has  led  Ashour-Abdalla  and  Kennel^10^ 
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to  derive  expressions  for  anomalous  velocity  diffusion  coeffi¬ 
cients.  Such  as  analysis  may  be  pertinent  in  (1)  describing 
the  fluctuation  level  of  the  magnetosphere  environment  through 
which  a  rocket  fired  beam  propagates,  or  (2)  suggesting 
how  a  steady  state  energetic  particle  distribution  created 
by  the  rocket  fired  beam  may  become  unstable. 

2.3.4  RF  Signals 

The  plasma  and  electron  cyclotron  frequencies 
are  expected  to  be  most  characteristic  of  the  rf  signals 
obtained  from  the  beam  plasma  interaction.  The  most  complete 
published  data  on  rf  signatures  of  the  beam  plasma  interaction 
in  the  ionosphere  are  due  to  Cartwright  and  Kellogg^^, 
who  analyzed  the  ECHO  I  experiment.  They  show  that  signals 
at  the  ambient  plasma  frequency  as  well  as  strong  bands 
of  radiation  at  about  twice  the  electronc  cyclotron  frequency 
were  emitted.  In  addition,  there  are  strong  sources  of 
radiation  below  the  electron  cyclotron  frequency. 

The  measurements  of  the  rf  signals  were  made  from 
a  nose  cone  antenna,  which  was  separated  from  the  main  rocket. 
The  antenna  was  therefore  immersed  in  the  ionospheric  plasma 
and  obtained  signals  from  the  beam  at  ever  increasing  distances 
during  the  flight.  The  need  for  a  nearby  antenna  for  rf 
detection  seems  to  be  a  necessary  requirement  for  many  rocket 
fired  beam  experiments.  Ground  based  rf  antennas  were 
generally  able  to  receive  signals  associated  with  only  a 
limited  number  of  the  radiation  sources  detected  by  the 
nose  cone  antennas. 

The  observed  plasma  frequency  radiation  was  computed 
by  Cartwright  and  Kellogg  as  incoherent  electrostatic  Cerenkov 
radiation,  and  the  value  obtained  was  in  agreement  with 
the  measurements.  However,  the  amplitude  of  the  signal 
unexpectedly  peaked  at  about  5  ms  into  the  16  ms  pulse 
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characteristic  of  ECHO  I.  In  addition,  the  peak  radiation 
occurred  for  those  injected  pitch  angles  just  upfield  from 
90°.  It  was  also  expected  that  the  radiation  would  drop 
off  as  the  inverse  of  the  square  root  of  the  distance  between 
the  nose  cone  antenna  and  the  beam.  The  amplitude,  however, 
did  not  vary  according  to  this  relation  as  the  nose  cone 
and  beam  separated  from  each  other. 

The  observations  of  whistler  modes  (f<f  )  also 

Cc 

led  to  surprising  results.  Near  apogee,  a  signal  extending 
down  only  to  a  frequency  of  about  f  /2  'v  500  kHz  was  obtained, 
while  at  lower  altitudes,  whistlers  were  observed  with  fre¬ 
quencies  down  to  50  kHz,  the  frequency  limit  of  the  receiver. 

In  addition,  the  radiation  amplitude  peaked  within  the  beam 
pulse  width  just  as  with  the  plasma  frequency  waves.  Again, 
the  maximum  amplitudes  occurred  for  a  pitch  angle  range 
just  up  field  from  90°. 

The  radiation  near  2  f  also  has  puzzling  features. 

C  V 

It  is  well  known  that  for  f  ^  >  f  ,  there  is  a  cold  plasma 

stop-band  above  fce  that  would  preclude  the  existence  of 

the  lower  harmonics  of  the  electron  cyclotron  frequency. 

Thus,  one  would  seem  forced  to  explain  the  2  fce  in  terms 

of  warm  plasma  theory.  Bernstein  modes  result  from  a  warm 

plasma  analysis  and  do  have  frequencies  near  the  electron 

cyclotron  harmonics.  However,  Cartwright  and  Kellogg  point 

out  that  the  Bernstein  mode  could  explain  the  relatively 

large  bandwidth  of  the  2  f  radiation  only  if  the  wavelength 

became  very  short  and  that  such  wavelengths  could  not  be 

observed  by  the  nose  cone  antennas.  Further  evidence  of 

the  second  electron  cyclotron  harmonic  has  been  reported 

( 12 ) 

by  Monson  and  Kellogg  who  carried  out  ground  based  rf 

antenna  observations  in  conjunction  with  the  ECHO  IV  experi¬ 
ment.  Detection  of  these  signals  on  the  ground  is  significant 
because  it  implies  an  efficient  conversion  process  through 
which  a  radio  wave  is  produced. 
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The  ARAKS  experiments  utilized  a  much  larger  beam 
current  than  in  the  case  of  ECHO  I.  Very  much  the  same 
kinds  of  signals  were  observed  by  the  nose  cone  antenna, 
except  that  the  predominant  electron  cyclotron  harmonies 
appeared  to  be  n=4  rather  than  n=2 .  One  important  difference 
in  the  experiments  was  that  the  signals  in  ARAKS  did  not 
peak  within  the  pulses  but  rather  were  continuous. 

Strong  ground  based  rf  signals  were  detected  from 

(13) 

the  ARAKS  and  ZARNITSA-2  experiments  (Cambou  et  al  ; 

( 14 ) 

Mishin  an  Ruzhin'  )  with  frequencies  in  the  tens  of 

MHz  range*  This  result  is  surprising  since  the  largest 

frequency  associated  with  the  plasma  is  the  upper  hybrid 

which  in  the  E  and  F  layers  takes  on  a  value  of  about 

3  MHz*  In  each  of  the  experiments,  the  signal  is  attributed 

to  the  production  of  a  high  density  plasma  in  the  vicinity 

of  the  rocket.  Creation  of  a  plasma  density  about  two 

orders  of  magnitude  larger  than  the  ambient  density  could 

explain  the  high  frequency  of  the  signals.  An  intensive 

analysis  of  the  physics  of  plasma  density  enhancement  has 

been  carried  out  by  Mishin  and  Ruzhin  with  respect  mainly 

to  the  ZARNITSA-2  expferiment .  According  to  Mishin  and  Ruzhin, 

the  trajectory  of  the  rocket  was  sufficiently  low  that  the 

condition  co  / cu  <1  was  valid  over  nearly  the  entire  flight, 
pe  ce 

These  authors  propose  for  this  case  that  the  cyclotron  reso¬ 
nance  branch  of  the  high  frequency  electron  mode  dispersion 
relation  is  ultimately  responsible  for  the  large  plasma 
density  enhancement.  They  argue  that  the  fluctuating  electric 
fields  heat  up  the  ambient  electrons  sufficiently  to  ionize 
the  neutral  particles  present  in  the  ionosphere.  This  mecha¬ 
nism,  whereby  plasma  electrons  are  created  by  ionization 
of  neutrals  from  ambient  electrons  heated  by  a  beam  plasma 
instability*  is  referred  to  as  the  beam-plasma  discharge 
(BPD).  The  magnitude  of  the  accelerating  turbulence  results 
from  saturation  of  the  linear  modes.  Mishin  and  Ruzhin 
have  estimate  a  steady  state-turbulence  level  by  appealing 
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to  the  arguments  of  strong  turbulence  theory  as  well  as 
considering  those  levels  appropriate  to  the  weak  turbulence 
(quasi-linear )  theory.  In  the  former  case,  the  theory  of 
weak  turbulence  breaks  down  because  of  a  modulational  insta¬ 
bility.  The  mechanics  of  the  stabilization  of  this  parametric 
mode  to  a  final  state  including  cavitons  is  not  well  understood 
at  this  time,  but  the  field  is  very  active.  It  has  yet 
to  be  determined,  conclusively,  whether  strong  turbulence 
is  responsible  for  naturally  occurring  rf  phenomena.  Of 
particule  importance  in  the  theory  is  that  a  critical  neutral 
number  density  exists,  above  which  more  particles  are  created 
through  ionization  processes  than  are  lost  by  diffusion. 

The  value  for  this  number  allows  for  an  experimental  check 
to  the  theory. 

Having  discussed  the  conditions  for  BPD ,  Mishin 
and  Ruzhin  finally  determine  the  conversion  efficiency  of 
the  electrostatic  mode  (basically  Langmuir  oscillations) 
to  a  radio  wave  that  can  be  received  at  a  ground  based  antenna. 
Two  processes  have  been  considered.  One  is  wave  fusion 
of  the  longitudinal  wave  with  an  ion  sound  mode  ( Z  +  s  -*■  t ) . 

The  other  suggested  process  is  the  fusion  of  the  Langmuir 
wave  with  still  another  longitudinal  wave,  the  Bernstein 
mode  ( l  +  l  -*■  t ) . 


The  HF  results  from  the  ZARNITSA-2  experiment 
as  well  as  the  theory  for  ti)pe<  wce  are  of  particular  inter¬ 
est  to  American  investigators  since  BPD  ground  experiments 

in  which  an  electron  gun  is  fired  into  a  gas  filled  chamber 

( 15  ) 

have  been  carried  out  by  Bernstein  et  al.  An  ambient 

magnetic  field  was  established  within  the  chamber  and  the 

pitch  angle  of  the  injected  beam  was  about  20°.  The  chamber 
-5  -6 

pressure  of  10  -  10  T  as  well  as  the  magnetic  field 

strength  of  just  over  1  G  are  indicative  of  conditions  in 
the  lower  ionosphere.  The  beam  energy  of  ^  1  keV  was  much 
lower  than  that  characteristic  of  rocket  fired  beam  experiments. 
Since  the  background  gas  was  not  pre-ionized,  an  ambient 
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backgroun  plasma  had  to  be  created  by  direct  ionization 

from  the  beam.  Thus,  for  low  beam  currents  (prior  to  BPD), 

the  condition  u>  <  was  satisfied.  The  experiments 

showed  that  at  a  critical  current,  the  electron  number  density 

increased  dramatically.  It  is  not  yet  clear  that  the  BPD 

chamber  experiments  have  been  adequately  explained.  Bernstein 
(15) 

et  al*  have  compared  the  experimental  results  with  a 

( 16  } 

theory  by  Lyachov  and  Manadadze^  ' ,  who  have  done  work 

on  the  ZARNITSA-2  experiment ,  and  have  concluded  that  the 

oredicted  critical  neutral  density  was  an  order  of  magnitude 

higher  than  required  in  the  experiment.  Rowland  and 
( 17 ) 

Papadopoulos  have  postulated  that  BPD  for  <  wce 

is  triggered  by  the  onset  of  the  Cerenkov  beam  mode  coupled 
with  the  slow  cold  electron  plasma  mode.  A  critical  ambient 
plasma  density  is  required  for  onset  of  the  mode.  Rowland 
and  Papadopoulos  claim  that  their  model  explains  the  experi¬ 
mental  scaling. 

BPD  appears  to  occur  when  beams  are  fired  within 

the  E  and  F  layers  where  u>  >  w  This  case  has  been 

pe  ce 

treated  in  less  detail  then  for  u>  <  <o  ,  and  no  ground 

pe  ce 

based  chamber  experiments  have  as  yet  been  carried  out  in 
this  regime. 

2.3.5  Turbulence  and  Anomalous  Beam  Transport 

Rocket  fired  electron  beam  experiments  have  given 
ample  evidence  of  the  presence  of  important  collective  effects. 
Some  of  the  observed  phenomena  which  cannot  be  explained 
on  the  basis  of  particle-particle  interaction  are:  enhanced 
local  backscatter  in  upward  injection;  diffuse  streak  effects 
in  downward  injection;  excessive  particle  flux  losses  in 
echoing  experiments;  and,  long  delay  times  in  the  return 
of  particles  that  have  mirrored. 
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effect  will  eventually  heat  the  plasma  so  that  the  mode 

would  continue  to  grow  when  kinetic  considerations  are  taken 

into  account.  It  has  long  been  felt  that  a  final  saturation 

level  could  be  formed  by  the  method  of  weak  turbulence, 

but  this  no  longer  may  be  the  case.  It  has  been  determined, 

for  example,  that  Langmuir  turbulence  cannot  be  described 

2 

by  weak  turbulence  theory  whenever  W/nT  >  (KX^)  ,  where 

W  is  the  turbulent  energy  density,  T  is  the  plasma  temperature, 

( 23 ) 

and  XD  is  the  Debye  length.  The  subsequent  behavior 

of  this  strong  or  spiky  turbulence  has  not  been  well 
established. 
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In  mirroring  (or  echo)  experiments  Where  the  beam 
propagates  within  the  magnetosphere,  it  is  important  to 
consider  the  turbulent  spectrum  caused  by  natural  effects. 

In  this  regard,  an  interesting  analysis  has  been  carried 
out  by  Ashour-Abdal la  and  Kennel  These  authors  assumed 

the  presence  of  a  stationary  spectrum  based  on  spatial  growth 
profiles  existing  in  the  inhomogeneous  magnetosphere.  The 
basic  linear  instability  considered  was  the  upper  hybrid 
electrostatic  electron  mode  being  driven  by  electron  cyclotron 
resonances  associated  with  a  hot  component  with  loss  cone 
anisotropy.  The  analysis  finally  allows  for  determination 
of  a  turbulent  spectrum  which  would  be  part  of  the  ambient 
background  through  which  a  rocket  fired  electron  beam  would 
propagate.  There  is  a  large  literature  on  magnetospheric 
structure,  and  this  is  an  area  of  importance  in  beam  propa¬ 
gation  studies.  A  recent  representative  work  in  this  field 

( 18 ) 

is  a  paper  by  Kan  and  Akasofu  ,  who  discuss  the  structure 

of  the  magnetosphere-ionosphere  required  to  support  an  auroral 

electric  field*  The  role  of  the  whistler  mode  in  magneto- 

spheric  turbulence  was  first  emphasized  in  an  important 

(19) 

paper  by  Kennel  and  Petschek  ,  and  their  notion  of  a 

flux  limit  on  stably  trapped  particles  in  the  magnetosphere 

( 20  21 ) 

has  recently  attracted  much  attention.  ’ 

As  discussed  above,  determination  of  the  turbulent 
spectrum  for  beam  plasma  interaction  must  proceed  mainly 
on  theoretical  grounds.  There  is  a  rather  good  justification 
for  examining  the  high  frequency  electrostatic  electron 
mode  driven  by  beam  Cerenkov  or  cyclotron  interaction  as 
being  response  for  self-generated  turbulence.  Analysis 
of  this  mode  for  cold  beam  has  been  examined  in  some  detail 
by  Papadopoulous  for  finite  systems  in  the  respective 

limits  of  large  magnetic  field  and  large  plasma  density, 
respectively.  On  the  other  hand,  the  mechanisms  responsible 
for  saturation  of  the  linear  mode  are  not  well  understood. 

The  beam  may  become  trapped  within  its  own  field,  but  this 
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effect  will  eventually  heat  the  plasma  so  that  the  mode 

would  continue  to  grow  when  kinetic  considerations  are  taken 

into  account.  It  has  long  been  felt  that  a  final  saturation 

level  could  be  formed  by  the  method  of  weak  turbulence, 

but  this  no  longer  may  be  the  case.  It  has  been  determined, 

for  example,  that  Langmuir  turbulence  cannot  be  described 

2 

by  weak  turbulence  theory  whenever  W/nT  >  (KA^)  ,  where 

W  is  the  turbulent  energy  density,  T  is  the  plasma  temperature 
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and  XD  is  the  Debye  length.  The  subsequent  behavior 

of  this  strong  or  spiky  turbulence  has  not  been  well 
established. 


Section  3 


RESULTS  AND  THEIR  INTERPRETATION  WITH  RESPECT  TO  ECHO-4  DATA 


3.1  TEST  RUNS  OF  MCBE 

We  have  carried  out  important  test  runs  of  MCBE, 
a  Monte-Carlo  algorithm  for  the  transport  of  multi-keV 
electrons  along  magnetic  field  lines.  The  code  has  been 
used  to  simulate  downward  rocket-fired  electron  beam 
scattering  and  deposition  in  the  ionosphere.  Using 
Gouldsmit-Saunderson  multiple  scattering  theory,  we  have 
determined  the  beam  width  in  the  region  of  energy  de¬ 
position.  Identification  of  these  results  with  beam  width 
measurements  of  artificial  aurora  shows  that  the  model 
accurately  computes  the  beam  width  resulting  from  particle- 
particle  effects.  We  believe  that  MCBE  gives  the  best 
available  predictive  capability  for  classical  beam  widths 
as  functions  of  the  requisite  beam  parameters — energy  and 
pitch  angle. 

Several  test  runs  were  carried  out.  Figure  2 

shows  the  energy  deposition  of  20  keV  electrons  injected 

downward  into  the  ionosphere  from  300  km  at  0°  pitch  angle 

as  a  function  of  altitude.  The  dotted  curve  represents 

24 

deposition  results  from  Berger  et  al .  In  both 
cases,  maximum  deposition  occurs  at  about  100  km  and  the 
results  are  in  good  quantitative  agreement.  A  similar 
plot  is  shown  in  Figure  3  where  deposition  is  into  a  N„ 
atmosphere,  and  the  energy  distribution  is  10  keV  gaussian. 
The  broken  curve  gives  results  from  B3C,  an  algorithm 
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ENERGY  DEPOSITION  (eV/cm) 


FIGURE  2.  Energy  Deposition  as  a  Function  of 
Altitude  for  Downward  Injection  of  20  keV  Electrons 
from  300  km  at  Normal  Incidence 


The  solid  histogram  is  the  MCBE  result,  while 
the  dotted  curve  gives  the  result  of  Berger, 
et  al£4 


ALTITUDE  (km) 


ENERGY  DEPOSITION  (eV/cm) 


FIGURE  3.  Energy  Deposition  as  a  Function  of 
Altitude  for  Downward  Injection  at  Normal  Incidence 


Energy  Distribution  is  10  keV  Gaussian.  The  solid 
histogram  is  -he  MCBE  result,  while  the  broken  curve 
is  the  result  from  the  B3C  algorithm. 


B 


developed  by  SAI  to  give  optical  emissions  from  electron 
deposition.  In  this  case,  MCBE  was  run  with  1000  histories 
and  computed  26%  baokscatter  across  the  injection  plane. 

The  corresponding  percentage  for  B3C  was  24%.  The  effect 
of  initial  pitch  angle  on  the  deposition  profile  is  shown 
in  Figure  4  for  20  keV  electrons.  The  light  solid  curve 
represents  an  isotropic  distribution  of  downward  directed 
pitch  angles.  Results  for  45°  pitch  angle  injection  are 
given  by  the  dotted  curve,  while  the  deposition  of  normally 
incident  electrons  are  shown  by  the  heavy  solid  curve. 
Deposition  profiles  for  normal  incidence  for  a  family  of 
electron  energies  are  presented  in  Figure  5.  Results  for 
energies  of  10,  20,  and  30  keV  are  plotted. 

Figure  6  presents  the  effective  beam  width  as  a 
function  of  altitude  for  20  keV  electrons  injected  normally 
(solid  curve)  and  isotropically  (dotted  curve).  The  inter¬ 
esting  point  here  is  that  electron  scattering  substantially 
spreads  out  normally  incident  electrons  so  that  at  the 
lower  altitudes  where  deposition  and  optical  emissions 
occur,  the  beam  spread  is  effectively  independent  of 
initial  pitch  angle  distribution.  Note  that  the  effective 
beam  radius  is  less  than  the  characteristic  quantity 
2|^|/wce  that  gives  a  measure  of  the  maximum  gyro-radius. 
This  result  will  be  discussed  further  in  the  next  section. 

3.2  INTERPRETATION  OF  ARTIFICIAL  AURORA 

EXPERIMENT  IN  ECH0-4  ROCKET  LAUNCH 

3.2.1  Comparison  with  MCBE  Results 

The  generation  of  artificial  aurora  in  ECHO  4 

p  c 

rocket  launch  (Hallinan  et  al.  °)  serves  as  a  relevant 
test  for  the  analyses  developed  in  the  AFGL  program.  A 
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ENERGY  DEPOSITION  (eV/cm) 


FIGURE  4*  Energy  Deposition  as  a  Function  of 
Altitude  for  Downward  Injection  of  20  keV  Electrons 

from  300  km 


The  light  solid  curve  represents  an  isotropic 
distribution  of  downward  directed  pitch  angles. 

The  dotted  curve  represents  45°,  while  the  heavy 
solid  curve  shows  the  results  for  normally  incident 
electrons . 


ENERGY  DEPOSITION  (eV/cm) 


FIGURE  5.  Energy  Deposition  as  a  Function  of 
Altitude  for  Downward  Injection  from 
300  km  of  Electrons  of  Various  Energies 


The  dotted  histogram  shows  results  for  10  keV  electrons 
the  heavy  solid  histogram  indicates  results  at  20  keV; 
the  light  solid  histogram  shows  results  for  30  keV 
electrons.  In  all  cases,  the  pitch  angle  is  0°. 
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FIGURE  6.  Effective  Beam  Width  as  a  Function  of 
Altitude  for  Downward  Injection  from 
300  km  of  20  keV  Electrons 


The  solid  histogram  shows  results  for  normal  incidence, 
while  the  dotted  histogram  gives  results  for  an 
isotropic  pitch  angle  distribution. 


capsule  description  of  the  ECHO  series  of  experiments  is 
given  in  the  previously  presented  Tables  1-5  and  in  Figure  1. 
The  experiment  was  particularly  interesting  because  the  radial 
extent  of  the  beam  could  be  resolved.  A  summary  of  the 
streak  width  results  is  given  in  Table  6. 

Optical  emissions  from  the  artificial  aurora 
emanated  from  the  region  of  maximum  energy  deposition, 
between  100  and  90  km.  As  shown  in  Table  6,  the  observed 
streak  widths  depended  strongly  on  the  altitude  at  which 
the  gun  was  fired.  At  rocket  apogee,  streak  diameters  well 
in  excess  of  100  m  were  found,  and  these  were  called  diffuse 
aurora.  All  other  examples  were  referred  to  as  normal 
streaks,  although  it  is  seen  that  there  is  a  considerable 

variation  in  the  streak  widths  at  the  other  firing  altitudes. 

2  5 

Hallinan  et  al.  interpreted  the  diffuse  aurora  as  being 

the  result  of  wave-particle  effects,  but  they  generally 

attributed  the  other  test  results  to  classical  effects.  We 

believe  that  their  interpretation  followed  from  the  results 

24 

of  Berger  et  al .  who  presented  beam  widths  in  excess  of 
100  m  due  solely  to  particle-particle  effects.  The  corre¬ 
sponding  results  from  MCBE,  as  shown  in  Figure  6,  give 
much  smaller  streak  widths  due  to  particle-particle  inter¬ 
actions.  The  width  of  the  beam  at  auroral  altitudes  is 
between  32  and  33  m,  in  good  agreement  with  the  ECHO  4  streak 
width  observed  for  a  firing  altitude  of  148  km.  Although 
the  energies,  as  well  as  the  pitch  angle  ranges,  for  this 
firing  are  not  the  same  as  in  the  MCBE  run,  the  signifi¬ 
cance  of  the  numerical  results  is  clear.  The  streak 
widths  due  to  particle-particle  interactions  are 
relatively  small,  and  significant  increases  in  the  widths 
must  be  due  to  other  effects,  such  as  wave-particle 
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TABLE  6.  SUMMARY  OF  ECHO  4  STREAK  WIDTH  OBSERVATIONS 


FIRING  ALTITUDE 
(km) 

ENERGY 

(keV) 

A 

r 

e 

Cm) 

STREAK  WIDTH 
(m) 

PITCH  ANGLE 
(°) 

148  (UPLEG) 

26 

11.8 

34 

70 

160 

36 

13.9 

44 

70 

210 

26 

11.8 

182 

70 

214 

38 

14.3 

128 

30-90 

181  (DOWNLEG) 

34 

13.5 

72 

70 

156 

33 

13.3 

59 

30-90  J 
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interaction.  We  therefore  conclude  that  wave-particle 
interaction  plays  a  role  in  firings  between  148  km  and 
apogee.  The  effects  are  continuous,  becoming  stronger 
as  the  altitude  increases. 

3.2.2  Wave-Particle  Effects  in  Artificial  Aurora 

Any  study  of  wave-particle  effects  involves  the 
accurate  specification  of  equilibrium  conditions.  In 
rocket  fired  e-beam  experiments,  we  have  found  that  one  of 
the  major  problems  is  beam  characterization.  In  most 
experiments,  the  electrons  are  injected  at  a  particular 
pitch  angle.  For  such  cases,  the  single  particle  tra¬ 
jectories  fall  on  a  helix  wound  about  a  cylinder  along 
the  magnetic  field  line  with  Larmor  radius  Vj_/o)  ,  where 
Vj_  is  the  electron  velocity  component  perpendicular  to  the 
magnetic  field  and  u>cp  is  the  electron  cyclotron  frequency 
(eB/mc  in  cgs  units).  A  pulse  of  electrons  would  occupy 
a  finite  length  sector  of  the  helix.  Since  the  pulse  width 
is  much  larger  than  the  time  for  an  electron  to  complete  a 
cyclotron  orbit,  the  axial  length  of  a  single  pulse  is  of 
order  thousands  of  time  longer  than  the  Larmor  radius. 

Space  charge  effects  can  change  this  picture  however.  If 
the  number  density  of  the  beam  at  the  gun  orifice  is  larger 
than  the  ambient  plasma  density,  the  beam  will  expand  at 
least  until  the  beam  density  reduces  to  the  ambient  density, 
at  which  point  the  charge  begins  to  be  neutralized.  This 
will  give  rise  to  velocity  dispersion.  For  example, 
following  Gendrin  ,  we  can  estimate: 


where  I  is  space  charge  limited  current.  For  experiments 
of  interest  where  E  =  25  keV,  I  =  131A.  For  ECHO  4 
currents  70  mA,  this  yields  a  2%  velocity  dispersion,  which 
for  pitch  angles  of  60°  is  sufficient  to  "fill  in"  the 
space  between  the  initial  helical  windings  in  about  1  mile. 
Of  course,  even  this  scenario  can  be  made  more  complex  by 
efforts  occuring  closer  to  the  rocket,  namely  instabilities 
associated  with  the  helical  tube  which  may  be  associated 
with  beam  plasma  discharge  (BPD).  At  any  rate,  it  seems 
reasonable  to  assume  that  the  beam  can  be  characterized  by 
reasonably  symmetric  shells  within  a  few  kilometers  of  the 
injection  point.  If  this  configuration  gives  rise  to  an 
instability,  then  the  electric  field  that  develops  opposes 
the  electron  motion  (i.e.,  reducing  the  beam  kinetic  energy). 
This  means  that  E  x  B  is  directed  outward  so  that  the 
electrons  will  drift  outward,  thus  increasing  the  radius 
of  the  beam. 

It  has  been  well  established  that  the  hollow  beam 
configuration  to  which  the  helix  has  evolved  is  unstable  to 
two  stream  modes.  In  the  simplest  two  stream  analysis,  the 
streaming  velocity  is  parallel  to  the  magnetic  field.  The 
hollow  beam  development,  however,  is  due  to  the  transverse 
anisotropy  of  the  electron  pulses.  This  configuration  can 
be  shown  to  produce  instability  analogous  to  the  usual  two 
stream  result.  The  nonlinear  saturation  of  the  cold 
electrostatic  wave  is  due  to  particle  trapping.  The  field 
amplitude  at  saturation  scales  as 


E? 

8t r 


aNbeb 


where  eb  is  the  beam  particle  energy  and  a  is  a  dimension¬ 
less  factor  between  about  0.1-0. 3.  The  magnitude  of  the 
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E  x  B  drift  of  resonant  particles  is  therefore 


V  =  =  vh 

Dm  w  b 

ce 


where  is  the  beam  particle  speed. 


We  now  compute  the  time  required  for  the  beam  to 
spread  out  to  a  radius  of  91  m,  the  largest  streak  width 
observed  in  the  ECHO  4  experiments.  Typical  beam  para¬ 
meters  for  the  experiments  are  electron  energies  of  25  keV 
and  pitch  angles  of  60°.  For  simplicity,  we  assume  that 
the  beam  is  solid  rather  than  hollow  and  has  an  initial 
radius  equal  to  the  Larmor  radius.  Since  the  Larmor  radius 
is  only  10m,  consideration  of  a  solid  cylinder  cannot 
change  the  result  significantly.  The  beam  plasma  frequency 


where  is  the  initial  beam  plasma  frequency  correspond¬ 

ing  to  the  cylinder  of  radius  rL.  The  equation  for  the 
expanding  radius  is  just 


dr 

dt 


Letting  a  =  0.1  and  B  =  0.46G,  we  find  that  the  beam  width 
evolves  to  91  m  radius  in  a  distance  of  17  km  from  the  gun 
injection  point. 

In  order  for  the  drift  to  be  sustained,  the 
system  must  remain  collisionless.  This  means  that  the 
higher  energy  electrons  may  not  undergo  collisions  with 


neutrals  during  the  wave  growth  period.  The  growth  rate  for 
many  modes  of  interest  is  of  order 

Y  %  u>p  =  1.6  •  106  sec-1 

O 

For  a  nighttime  density,  nD  =  2  *  10°  at  altitudes  between  100 
and  200  km,  the  growth  rate  for  a  200m  diameter  beam  with 
the  aforementioned  parameters  is 


y 


1.33 


10' 


-1 

sec 


These  growth  rates  are  much  higher  than  the  primary  electron 
collision  frequencies.  In  fact,  25  keV  electrons  experience 

5 

a  collision  frequency  with  neutrals  of  only  about  1.8  •  10  sec 
at  an  altitude  as  low  as  130  km.  Hence,  for  altitudes  above 
this,  the  beam  plasma  system  is  effectively  collisionless. 

The  question  arises  as  to  why  even  greater  beam  expansions 
than  indicated  by  the  maximum  observed  diameters  of  ,v200m 
do  not  occur.  The  reason  may  be  that  the  type  of  mode 
changes.  Note  that  the  plasma  density  for  which  u>p  =  u>ce 
is  just 


n  =  2.06  •  104cm  3 

P 

For  nighttime  electron  densities  associated  with  artificial 

aurora  experiments,  the  plasma  density  may  be  smaller  than 

this  value  at  altitudes  of  200  km  and  below.  In  this 

region,  the  excited  modes  would  be  subject  to  the  condition 

u  /u)  >>1.  This  change  in  the  mode  topology  may  be  effec- 
ce  p 

tive  in  terminating  the  large  amplitude  fluctuation  spect¬ 
rum.  On  the  other  hand,  these  modes  may  still  be  respon¬ 
sible  for  streak  widths  greater  than  those  expected  from 
particle-particle  interactions  occurring  at  the  firing 
altitudes  under  200  km. 
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The  above  analysis  suggests  that,  diffuse  aurora 
(where  beam  widths  are  much  greater  than  100m)  will  occur 
whenever  the  plasma  frequency  is  larger  than  the  electron 
cyclotron  frequency.  Our  interpretation  of  why  beam  spread 
may  decrease  considerably  for  beam  propagation  under  200  km 


as  well  as  the  computed  spread  ratio  in  regions  where  oj  /wce  >  1 
are  consistent  with  experiment  where  diffuse  aurora  were 


obtained  for  beam  injection  altitudes  of  210  km  and  214  km. 


It  is  interesting  to  note  that  beam  spread  is  insensitive 


to  the  beam  number  density  (i.e.,  the  beam  current  for  given 
energy  and  pitch  angle).  For  a  given  distance  along  the 


field  line  over  which  the  instability  develops,  the  beam 
spread  varies  only  as  n^.  Such  insensitivity  of  the  width 
of  an  artificial  aurora  to  beam  current  can  easily  be 
checked  by  experiment. 


Other  factors  need  be  considered  before  the 
above  theory  can  indeed  be  verified.  First,  the  evolution 
of  the  helical  beam  to  a  neutral  hollow  cylinder  and  then, 
perhaps  a  solid  cylinder  needs  to  be  investigated  experi¬ 
mental  1 y,  relat ivel y  close  to  the  electron  gun.  Second, 
the  temperature  of  the  beam  should  be  determined.  If  the 
beam  heats  up  rapidly,  the  instability  will  shift  from  the 
hydrodynamic  "cold"  mode  to  a  kinetic  "hot"  mode.  Such 
modes  not  only  have  smaller  growth  rates  in  the  linear 
regime,  but  may  be  susceptible  to  caviton  formation  in 
the  nonlinear  regime.  It  has  been  suggested  that  cavitons 
could  accelerate  ambient  electrons  to  high  energies.  It 
would  bo  of  interest  to  measure  the  particle  spectrum  to 
determine  whether  this  effect  occurs.  A  third  factor  is 
that  the  changes  in  mode  structure  can  influence  the  con¬ 
vective  nature  of  the  instability,  so  that  whereas  a  cold 
plasma  mode  may  remain  effectively  stationary  in  radial 
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position,  a  warm  mode  may  convect  quickly  out  ol'  the 
generation  region  before  the  fluctuation  level  becomes 
large  enough  to  influence  particle  orbits.  The  problem 
of  mode  convection  within  warm  finite  beams  is  yet  to  be 
addressed . 


Section  4 


SUGGESTIONS  FOR  ROCKET  FIRED  BEAM  EXPERIMENTS 

AT  AFGL 

4.1  ION  BEAM  PROPAGATION 

One  of  the  interesting  results  from  the  spacecraft 

charging  experiment  of  Cohen  was  the  fact  that  the  return 

current  very  closely  equalled  the  emitted  8  pA  ion  current, 

irrespective  of  the  rocket  altitude.  At  the  high  altitudes, 

where  the  ambient  ion  density  is  high,  this  finding  is  not 

surprising  since  the  return  current  can  be  obtained  from  the  large 

reservoir  of  electrons  traveling  no  faster  than  the  thermal 

speed.  At  an  altitude  of  150  km,  however,  where  the  ambient 

3 

density  is  only  10  ,  this  argument  would  yield  a  return  current 
of  less  than  0.1  pA.  Thus,  the  soruce  of  the  measured  return 
current  at  low  altitudes  is  not  readily  apparent.  It  has 
been  often  suggested  that  the  fluctuating  electric  fields 
arising  from  plasma  instability  may  increase  the  number  density 
by  accelerating  the  secondary  electrons  to  energies  above 
the  threshold  for  ionization  of  the  background  neutral  particles. 
The  accompanying  enhancement  of  the  background  electron 
number  density  could  be  sufficient  to  produce  the  required 
return  current.  We  propose  that  the  r.f.  signatures  of  insta¬ 
bilities  caused  by  ion  beams  should  be  monitored. 

The  use  of  ion  beams  in  rocket  fired  beam  experi¬ 
ments  should  allow  for  the  excitation  of  VLF  electrostatic 
modes.  According  to  theory,  the  hydrodynamic  magnetized 
ion-ion  instability  should  be  present  for  the  low  velocity 
xenon  ion  guns  that  have  been  used  by  Cohen  on  a  sounding 


rocket  and  on  the  SCATHA  satellite.  The  mode  propagates 
almost  perpendicular  to  the  magnetic  field,  while  it  is  driven 
by  a  beam  directed  at  least  partly  across  the  magnetic  field 
line . 


The  major  condition  for  driving  the  mode  is  that 
the  beam  particle  speed  be  less  than  a  quantity  of  order 
the  Alfven  speed: 


V 


A 


B 

V4irnMi 


For  ionospheric  altitudes  where  0+  is  the  dominant  ion  species 
7 

=  4*10  cm/sec.  It  is  relatively  easy  to  create  an  ion 
beam  with  small  enough  particle  velocities.  Since  the  xenon 
atom  is  so  massive,  the  velocity  of  1  keV  beam  particles 
is  only  3.8*10°  cm/sec,  well  below  the  Alfven  speed  threshold. 
In  fact,  the  greater  challenge  is  to  install  on  a  rocket 
an  ion  gun  where  beam  particle  speeds  can  be  made  to  exceed 
the  Alfven  speed,  so  that  the  ion-ion  instability  can  be 
avoided. 


The  mode  frequency  will  be  near  the  lower  hybrid 

frequency,  which,  for  ionospheric  parameters  of  interest, 

is  about  equal  to  the  ambient  ion  plasma  frequency.  In  the 

region  where  0+  is  dominant,  and  where  we  take  an  average 
5  -3 

density  of  10  cm  ,  we  find  an  expected  mode  frequency 

5 

w  =:  10  rad/sec  or  f  =  17  kHz,  i.e.,  within  the  VLF  range. 

The  instability  growth  rate  is  also  of  order  of  the  lower 
hybrid  frequency  and  is  much  larger  than  typical  ion  collision 
frequencies.  Thus,  collisional  damping  effects  can  be  ignored 
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The  above  mode  is  magnetic  in  the  sense  that  the 
wave  number  must  be  much  smaller  than  the  inverse  electron 
Larmor  radius.  On  the  other  hand,  the  mode  is  stabilized 
when  the  wave  number  is  so  small  that  an  electromagnetic  term 
in  the  dispersion  relations  takes  on  a  large  value.  The 
first  condition  implies  that  the  wavelength  cannot  be  smaller 
than  about  10  cm  (where  an  ambient  temperature  of  0.1  eV 
has  been  assumed).  The  second  condition  can  be  shown  to 
limit  the  wavelength  to  tens  of  meters.  For  the  xenon  beam 
conditions  quoted  above,  the  maximum  growth  rate  for  the 
linear  instability  will  occur  for  a  wavelength  of  about  1  m. 

The  instability  can  be  expected  to  most  affect 

the  beam  properties  when  the  beam  density  is  comparable  to 

the  ambient  density.  If  the  beam  density  at  the  nozzle  is 

larger  than  the  ambient  density,  then  it  is  generally  accepted 

that  the  beam  will  spread  out  until  a  stationary  configuration 

is  achieved.  The  density  is  expected  to  decrease  to  near 

the  ambient  level  at  which  point  the  beam  would  be  charge 

neutralized.  Therefore,  the  resulting  beam  plasma  coupling 

would  be  strong.  The  level  of  ion  current  needed  to  trigger 

the  instability  would  be  that  value  for  which  the  beam  density 

just  equalled  the  ambient  density  at  the  gun  nozzle.  The 
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current  corresponding  to  a  beam  density  of  10  cm  entering 
a  5  cm  diameter  nozzle  is  just  over  a  pA.  Increase  of  the 
current  enhances  the  volume  over  which  the  interaction  would 
take  place  and,  therefore,  produces  a  stronger  signal.  It 
has  been  suggested,  however,  that  a  beam  that  is  initially 
non-neutral  may  attain  a  steady  unneutralized  configuration. 
For  the  case  of  electron  beam  experiments,  it  is  still  not 
clear  whether  the  beams  are  actually  neutral.  The  absence 
of  charge  neutrality  in  an  ion  beam  experiment  could  prevent 
the  strong  beam  plasma  coupling  we  have  been  discussing. 
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DETECTION  OF  ELECTROMAGNETIC  SIGNALS  ON 
NOSE  CONE  ANTENNA 


The  use  of  nose  cone  antennas  has  been  an  important 
feature  of  the  ECHO  experiments  in  the  reception  of  r.f. 
signals.  The  nose  cones  are  ejected  from  the  rocket  and 
have  been  used  to  monitor  signals  at  distances  of  a  few  meters 
to  about  1  kilometer  from  the  beam.  Since  it  is  generally 
supposed  that  electrostatic  signals  are  most  strongly  excited 
in  the  beam  plasma  interaction,  the  r.f.  antennas  in  the 
ECHO  series  have  been  capable  of  measuring  only  electric 
fields.  It  is  of  interest,  however,  to  monitor  the  presence 
of  electromagnetic  waves  with  the  nose  cone  antenna.  It 
should  be  possible  to  measure  the  fluctuating  magnetic  field 
with  two  small  lightweight  loop  antennas  oriented  at  right 
angles  to  each  other.  Not  only  will  this  diagnostic  allow 
for  the  pickup  of  direct  electromagnetic  signals  from  the 
region  of  beam  plasma  interaction,  but,  perhaps  more  importantly, 
it  will  allow  for  the  detection  of  electrostatic  modes  that 
are  converted  to  electromagnetic  signals.  The  conversion 
mechanism  of  the  plasma  mode  to  a  vacuum  signal  is  well  known. 

In  addition,  however,  the  ECHO  experiments  have  picked  up 
signals  at  twice  the  electron  gyrof requency  on  the  ground, 
as  well  as  on  electrostatic  nose  cone  antennas.  This  mode 
has  generally  been  identified  as  the  Bernstein  mode  which 
propagates  nearly  perpendicular  to  the  ambient  magnetic  field. 

At  present,  the  conversion  mechanism  is  not  known.  The  ability 
to  detect  an  electromagnetic  signal  at  twice  the  gyrof requency 
(—3  MHz)  would  permit  the  determination  of  whether  or  not 
conversion  occurs  in  the  vicinity  of  the  rocket.  Even  a 
negative  result  here  would  be  of  interest,  because  it  would 
prove  that  the  conversion  took  place  away  from  the  rocket 
at  some  well-defined  position  in  the  ionosphere.  There  are, 
for  example,  theories  suggesting  that  conversion  takes  place 
because  of  ambient  density  gradients. 
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We  believe  that  Langmuir  probes  should  be  attached 
to  the  nose  cone  to  monitor  the  fluctuation  level  of  those 
modes  that  propagate  away  from  the  interaction  region.  Since 
the  characteristics  of  modes  propagating  within  a  stable 
plasma  away  from  the  source  region  are  fairly  well-understood 
in  terms  of  group  velocities  and  damping  rates,  it  should 
be  possible  to  reconstruct  the  fluctuation  data  to  evaluate 
conditions  within  the  interaction  zone.  Ed  Szuzczesicz  wants 
to  attempt  fluctuating  measurements  on  the  shuttle.  He  informs 
us  that  the  power  supply  required  for  the  modulated  signal 
on  the  probes  weigh  only  five  pounds.  It  thus  may  be  possible 
to  incorporate  such  a  diagnostic  on  a  rocket  nose  cone. 

4.3  CURRENT  NEUTRALIZATION  AND  ROCKET 

CHARGING  EXPERIMENTS 

As  discussed  before,  one  of  the  important  aspects 
of  beam  ejection  from  rockets  is  the  avoidance  of  potential 
build-up.  It  is  easy  to  show  that  in  the  absence  of  return 
current,  even  moderate  current  outflow  will  rapidly  increase 
the  vehicle  potential  so  as  to  shut  down  the  current.  Under 
many  conditions,  the  return  current  thought  to  be  available 
from  the  ambient  plasma  is  not  suffici  t  to  produce  the 
required  neutralization.  Interestingly,  however,  neutralization 
seems  to  occur  under  most  conditions  of  beam  injection.  The 
myriad  processes  that  must  cause  neutralization  are  not  under¬ 
stood,  and  we  suggest  a  program  to  examine  this  question. 

The  key  aspect  of  this  experimental  effort  would  be  to  modulate 
the  magnitude  of  the  beam  current.  As  the  current  increases, 
it  is  expected  that  the  following  sequences  of  neutralization 
processes  will  occur: 

•  flow  of  thermal  ambient  electrons 
toward  vehicle; 

•  potential  enhanced  flow  of  ambient  electrons; 

•  ionization  of  neutrals  to  enhance  electron 
number  reservoir  about  vehicle  (Beam  Plasma 
Discharge) . 
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With  currents  available  of  order  1A,  it  is  expected 
that  the  threshold  of  the  beam  plasma  discharge  (BPD)  will 
be  achieved.  Here  it  is  of  particular  interest  to  test  the 
ionospheric  scaling  and  compare  this  with  the  scaling  achieved 
at  the  Johnson  Space  Center.  The  results  of  the  laboratory 
experiments  suggest  the  BPD  is  not  very  sensitive  to  pitch 
angle,  and  this  is  another  result  that  this  experiment  could 
verify. 
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Determination  of  Artificial  Auroral  Streak  Widths  Using 
Guiding  Center  Diffusion  Theory 
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The  theory  of  two-dimensional  random  walk  is  applied  to  the  problem  of  guiding  center  diffusion  of 
kilovolt  electrons  caused  by  particle-particle  interactions.  The  mean  radius  of  the  guiding  center  distribu¬ 
tion  is  shown  to  depend,  in  a  simple  way,  on  the  scattering  and  energy  loss  properties  of  the  interactions. 

A  formula  for  the  width  of  an  artificial  auroral  streak  is  derived  which  is  consistent  with  observations  and 
Monte  Carlo  results. 


Introduction 

This  paper  addresses  the  problem  of  the  spreading  of  an 
electron  beam  under  the  influence  of  particle- particle  inter¬ 
actions  in  the  presence  of  a  magnetic  field.  The  specific  appli¬ 
cation  is  to  rocket-fired  electron  (e)  beams  in  the  ionosphere. 
The  basic  quantity  to  be  derived  is  the  average  radius  of  en¬ 
ergy  deposition.  The  work  is  analytic  and  considers  such  key 
parameters  as  the  Goudsmit-Saunderson  multiple  scattering 
distribution,  a  single  scattering  distribution,  and  the  probabil¬ 
ity  distribution  for  random  walking  in  two  dimensions.  The 
analytic  model  was  developed  to  provide  physical  insight  into 
results  being  generated  by  a  recently  developed  Monte  Carlo 
code  which  describes  the  transport  of  kilovolt  beam  electrons 
in  the  converging  geomagnetic  field  at  ionospheric  altitudes. 
A  paper  is  now  being  written  describing  the  capabilities  of  this 
code  and  results,  which  are  directly  comparable  with  experi¬ 
ment. 

In  addition  to  the  desire  for  physical  insight,  the  work  being 
reported  was  also  motivated  by  discrepancies  between  pre¬ 
dicted  beam  spreading  using  the  above  Monte  Carlo  code  and 
that  of  Berger  el  al.  [1970].  We  do  not  find  the  increase  in 
spreading  with  increasing  beam  energy  as  reported  by  them, 
although  altitude  profiles  of  energy  deposition  are  the  same. 
As  an  example,  we  predict  -]  as  much  spreading  at  a  starting 
beam  electron  energy  of  20  keV. 

The  spreading  of  e  beams  is  a  subject  of  considerable  inter¬ 
est  following  the  many  rocket  e  beam  experiments  over  the 
past  decade  (see  review  by  Winckler  [1980])  and  planned  ex¬ 
periments  from  both  rockets  and  orbital  vehicles.  Much  of  this 
interest  focuses  on  the  importance  of  wave-particle  effects  in 
spreading  beams  beyond  their  widths  produced  by  particle- 
particle  effects  alone,  i.e.,  by  multiple  scattering  and  energy 
loss  due  to  the  neutral  panicle  background.  There  is  evidence 
from  the  past  experiments  that  wave-particle  effects  can,  in 
fact,  be  disruptive  to  the  beam.  Perhaps  the  most  obvious  case 
comes  from  ECHO  IV  where  some  optically  measured  streaks 
are  diffuse,  while  others  appear  sharp  [ Halltnan  et  ai,  1978], 
As  much  as  a  factor  of  5  difference  in  widths  is  observed  for 
streaks  in  these  categories. 


Analytic  Model  and  Discussion 

The  applied  analytic  model  is  based  on  the  following  phys¬ 
ical  description.  An  electron  begins  with  energy  £„  directed 
parallel  to  the  magnetic  field.  Its  distribution  in  pitch  angle  af¬ 
ter  N  scatterings  is  given  by  the  Goudsmit-Saunderson  mul¬ 
tiple  scattering  distribution.  Energy  loss  is  assumed  to  be  con¬ 
tinuous,  given  by  the  Bethe  stopping  power  formula.  Finally, 
the  distribution  of  angles  for  any  single  scattering  event  is 
given  by  a  screened  Rutherford  formula.  The  model  may  be 
conveniently  divided  into  two  parts.  One  part  shows  that  the 
radial  aspect  of  the  electron  transport  can  be  represented  by 
two-dimensional  random  walking  of  the  guiding  center  (GC). 
The  key  parameter  to  come  from  this  part  is  /„  the  average 
GC  displacement  in  the  x-y  plane  (perpendicular  to  B)  caused 
by  the  Ah  collision.  The  other  part  of  the  model  applies  the 
concept  of  random  walking  in  two  dimensions  giving  the  GC 
radial  distribution  as  a  function  of  N,  the  number  of  scatter¬ 
ings.  From  this  comes  the  quantity  sought  in  this  work,  the  av¬ 
erage  radius  of  the  energy  deposition.  We  choose  to  present 
this  part  of  the  model  first,  which  leads  us  immediately  to  the 
results.  The  other  part,  some  of  it  being  of  a  geometric  nature, 
is  presented  in  Appendix  A. 

We  begin  with  the  guiding  center  located  at  r  -  0.  Upon  the 
Ah  collision  it  makes  an  average  displacement  I,  in  the  x-y 
plane  random  in  direction.  To  apply  the  GC  distribution 
function  below,  I,  must  be  constant.  Since  we  are  interested  in 
an  approximate  upper  bound  to  the  spreading  caused  by  par¬ 
ticle-particle  effects,  ln  which  is  a  slowly  varying  function  of  i, 
is  replaced  by  /„  an  upper  bound  to  /,  (see  Appendix  A).  The 
needed  probability  distribution  can  now  be  expressed  by 


W^r)  -  exp  - 


with  normalization  given  by 


f~  W*{r)2ardr-l 


1  Present  address  is  Bell  Laboratories,  HolmdeL  New  Jersey  07733. 
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Equation  (2)  is  the  two-dimensional  random  walk  distribution 
applicable  to  N  »  I  and  follows  from  the  three-dimensional 
expression  given  by  Chandrasekhar  [1943].  We  see  that  the 
distribution  becomes  flatter  and  wider  as  the  number  of  scat¬ 
terings  increases.  The  average  radius  of  the  particle  distribu- 


4763  The  U.S.  Government  H  authorised  to  reproduce  and  sail  this  report. 
A- 2  Permission  for  further  reproduction  by  others  must  be  obtained  rrom 
the  copyrlfht  owner. 


Paper  number  1 A0 1 84. 
0I48-O227/8I/0O1A-OI84SOI.QO 


4764 


Lin  and  Strickland:  Guiding  Center  Diffusion  Theory 


lion  after  N  scatterings,  <rN),  is 

<r„>«(»rlV),/1  y  (3) 

In  our  model,  energy  loss  is  continuous  along  the  electron 
trajectory,  which  provides  a  sound  description  for  energies  of 
interest,  i.e.,  keV  energies.  Thus,  between  scatterings  the  elec¬ 
tron  deposits  energy  uniformly  on  a  radius  r,  (electron  gyro- 
radius)  centered  about  the  given  guiding  center.  Let  us  assume 
that  this  deposition  takes  place  at  the  guiding  center.  Then  the 
distribution  for  energy  deposition  is  given  by 
K 

D(r)  -  2  A£,  W^r)  eV/unit  area  (4) 

(*■1 


where  N0  is  the  number  of  scatterings  corresponding  to  total 
energy  loss.  The  average  radius  of  energy  deposition  (r)D  is 

<r>»«  !  +  ■£<.  S  (5) 

t  /-i 


where  the  sum  comes  from  use  of  (4)  and  W \  is  given  by  (1). 
The  unit  length  is  taken  to  be  the  gyroradius  uju,  where 
and  <o  are  the  initial  speed  and  gyrofrequency.  H'  is  not  accu¬ 
rate  for  small  i,  and,  in  turn,  the  first  several  summed  terms  in 
(3)  are  not  accurate.  This  will  not  produce  significant  error, 
however,  because  of  the  large  size  of  N0.  The  unit  term  above 
accounts  for  deposition  taking  place  a  gyroradius  from  the 
guiding  center. 

We  can  cast  <  rD )  into  a  more  usable  form  noting  that  A£,  is 
slowly  varying  and  can  thus  be  removed  from  the  summation 
in  its  averaged  form  (see  Appendix  B).  We  now  have 


(r)n  a  1  + 


yfrr  KE 
2  £„ 


(6) 


Since  the  number  of  scatterings  N  -  £„/5E  is  large,  we  may 
take  the  continuum  limit,  giving 


</•)„»  1  + 


isX'1 

2  “£„  3|azj 


•  +  j 


(7) 


Inserting  the  expressions  for  5E,  and  from  the  appen¬ 
dices,  we  arrive  at  the  desired  form: 


(rlo*  1  +  J 


nZV'2 

In  EJE,  -  2 

2 

In  EJl  -  1 

(8) 


It  is  worth  noting  the  simple  physics  revealed  by  this  equa¬ 
tion.  The  numerator  inside  the  last  square  root  represents 
scattering,  while  the  denominator  represents  energy  loss.  If  ei¬ 
ther  is  set  to  zero,  which  corresponds  to  turning  off  that  mech¬ 
anism,  <r)c  gives  the  expected  result.  The  actual  value  of  the 
square  root  quantity  will  be  near  unity  for  our  applications, 
and  thus  the  average  radius  for  energy  deposition  is  about  two 
initial  gyroradii.  This  is  shown  in  Figure  1  along  with  avail¬ 
able  observations  for  rocket  e  beam  experiments.  We  note  that 
our  simple  form  for  <r)D  provides  reasonable  agreement  for 
beam  injection  at  low  attitudes.  Expression  (8)  and  our  Monte 
Carlo  results  are  independent  of  altitude,  in  contrast  to  the 
data.  The  spreading  observed  for  higher  injection  altitudes  is 


presumably  caused  by  a  combination  of  particle-particle  and 
wave-particle  effects. 

In  summary,  we  have  shown  that  the  lateral  spreading  of  an 
e  beam  by  particle-particle  interactions  can  be  formulated  as  a 
two-dimensional  random  walk  problem  involving  the  electron 
guiding  center.  Applying  the  concept  of  random  walking,  a 
simple  expression  is  derived  for  the  average  radius  of  energy 
deposition  for  an  initially  monoenergetic,  field-aligned  e 
beam.  This  expression  should  overestimate  rather  than  under¬ 
estimate  the  beam  spreading,  although  not  by  much.  The 
overestimate  is  based  on  the  use  of  /„  rather  than  I,  as  the 
radial  displacement  per  collision.  Subsequent  approximations 
that  were  made  should  not  bias  the  results  in  a  similar  fashion. 
We  find  that  the  average  radius  of  energy  deposition  is  only 
about  two  initial  gyroradii  and,  furthermore,  is  only  weakly 
dependent  on  energy.  The  observed  dependence  in  spreading 
is  in  conflict  with  the  Monte  Carlo  results  of  Berger  el  al. 
[  1970]  but  is  in  good  agreement  with  our  own  Monte  Carlo  re¬ 
sults.  We  further  observe  generally  good  agreement  with  the  e 
beam  spreading  measurements  of  Hallinan  et  al.  [1978]  for 
beam  injection  below  ~  160  km  and  with  the  Monte  Carlo 
models  of  McEntire  [1973]  and  Halbleib  [1974]. 

Appendix  A:  Derivation  of 
Radial  Displacement  /, 

The  concept  of  guiding  center  diffusion  is  introduced  to¬ 
gether  with  a  derivation  of  the  radial  displacement  of  the 
guiding  center  per  scattering.  Use  will  be  made  of  both  mul¬ 
tiple  scattering  and  single  scattering  concepts.  Consider  an 
electron  starting  with  energy  £„  moving  downward  along  a 
geomagnetic  field  line.  As  it  proceeds,  it  loses  energy  and  un¬ 
dergoes  scatterings,  thereby  providing  a  perpendicular  com¬ 
ponent  to  its  velocity  and,  in  turn,  precession  about  some 
fixed  guiding  center  between  scatterings.  Let  »,  be  the  velocity 
just  as  the  electron  begins  its  Ah  elastic  collision.  In  the  x  -  y 
plane  the  electron  is  positioned  at  (x„  y),  while  its  guiding 
center  is  at  (x,‘,  y,1).  These  coordinates  are  related  by 


x,‘  -  (Al) 

W 

and 

<A2) 

0 

where  u  is  the  gyrofrequency. 

The  electron  now  undergoes  its  Ah  collision.  While  its  posi¬ 
tion  does  not  immediately  change,  its  velocity  vector  does, 
and  in  turn,  so  does  the  position  of  the  guiding  center.  We 
wish  to  know  the  relationship  between  the  changes  in  these 
quantities.  Let  (u„  6, 4>)  be  the  spherical  coordinates  of  v„  in  a 
fixed  frame  with  field  B  along  z.  Further,  let  6  and  ^  be  the 
pitch  and  azimuthal  scattering  angles  with  respect  to  v  Then 
the  change  in  v*  in  terms  of  its  components  projected  onto  the 
x  -  y  plane,  can  be  shown  to  be 

do,'  ”  bJ-cos  0cos$sin0cos4  +  sin*sintfsind 

+  sin  0  cos  •  (cos  8  -  1)]  (A3) 
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Fig.  1.  Mean  radius  of  energy  deposition  as  a  function  of  height 
for  several  streaks  of  different  energies.  Numbers  above  data  points 
arc  beam  energies  in  keV. 


where  the  averaged  Legendre  polynomial  is 

(P1  cos  „  /^Xcos  gy/eos  8)  d  cos  8 
“**»/  ///cos  0)  decs  8 

For  the  single  scattering  distribution  /y  we  use 


//cos *)-/(£„ cos 9)- 


1 


(I  -  cos  8  +  2%)* 


(A9) 


(A  10) 


with  screening  parameter  t)y  given  by 


Expression  (A  10)  with  this  form  of  comes  from  use  of  a 
Thomas-Fermi  model  with  exponential  screening  by  the 
bound  electrons.  We  equate  E,  with  the  constant  terms  in 
(AM),  giving 

U,  =  E./E,  (A12) 

E,  appeared  explicitly  in  the  key  quantity  derived  in  this  pa¬ 
per,  (r)o,  and  thus  for  the  convenience  of  the  reader  we  note 
its  form  as  a  function  of  Z: 


a 


Z1'5 

0.885 


me' 

h1 


8£, 


(All) 


and 


E,  -  4.3 Z1'3  eV  (A  13) 


Ad,'  —  vH — cos  0  sin  sm  9  cos  <p  -  cos  9  sin  0  sin  <p 

Returning  to  p,  and  /,,  their  normalizations  are 

+  sin  0  sin  9  (cos  8  -  1)] 

Noting  from  (Al)  and  (A2)  that 

(A4) 

and 

J  pHca s  0)  2ir  d  cos  0  =»  1  (A14) 

Ad  ' 

A*/- - 1- 

(A5) 

J  IHcos0)2irdcos8~  1  (A15) 

and 

The  probability  density  function  for  the  Ah  step  to  have 
step  size  can  now  be  given.  It  is 

A  <  AW*‘ 

(A6) 

r  Hr)- 

2  /pXcos  0)/Xcos  8)  8  (r2  -  Ar,2)  dtp  d  cos  8  d  cos  8 

(A  16) 

we  are  now  able  to  relate  the  GC  radial  displacement  dr,  to 
the  change  in  velocity.  In  squared  form,  it  is 

Ar,1  -  ~  [cos1 9  sin2  9  cos2  <p 

+  sin2  8  sin2  <J>  +  sin2  8  (cos  8  -  l)2 
—  |  sin  20  sin  8  cos  <p  (cos  8  -  I))  (A7) 

We  note  that  Ar,2  is  independent  of  9,  which  introduces  the 
randomness  in  radial  motion  of  the  guiding  center  that  we  are 

seeking. 

To  proceed,  we  must  introduce  the  probabilistic  nature  of 
the  problem,  since  the  angles  in  (A7)  can  only  be  related  to 
the  Ah  scattering  through  probability  distribution  functions. 
Specifically,  a  multiple  scattering  distribution  function  gov¬ 
erns  the  behavior  of  8,  while  a  single  scattering  distribution 
Auction  governs  that  of  8.  A  uniform  distribution  is  assigned 
to  <p  due  to  the  azimuthal  symmetry  of  the  single  scattering 
process.  Let  pHcos  0)  and  /,( cos  0)  be  the  multiple  and  single 
scattering  distributions.  We  use  for  p,  a  generalization  of  the 
Goudsmtt  and  Sounder  son  (1940]  formula,  which  is 

a  m 

PH<xx0)-j-  X  (2f+  »)  n  <*Xcosl)>/Xcoae)  *>  I 

(A8) 


where  the  delta  function  provides  the  necessary  restrictions  on 
integrations  limits  such  that  any  given  combination  of  angular 
values  gives  only  r.  The  normalization  on  t  is 

2r  Jr,(r)r  dr  —  1  (A17) 

In  the  main  text  we  introduced  the  average  guiding  center 
displacement  for  the  Ah  collision  /,.  To  apply  the  given  two-di¬ 
mensional  random  walk  model  l,  was  replaced  by  its  value 
at  the  starting  energy.  There  we  noted  that  /„  provided  an  up¬ 
per  bound  to  I,.  We  now  give  the  form  of  /,  (in  squared  form) 
from  which  it  will  be  clear  that  /„  is  its  maximum  value.  It  is 
the  average  value  of  Ar2  and  is  thus  given  by 

I,1  m  2ir  JVr/frlr  dr 

-  2ir  I&r,1  pH  cos  9)/^cos  8)dcos0d<x>s8dp  (A18) 
In  gyroradius  units  referenced  to  E^vjui)  this  becomes 


f  (sin2tf),ii-tggglg>i 

Eo  2 


(A  19) 


where  the  ((cos  0  -  I)2),  terra  in  Ar,2  has  been  dropped.  Its 
contribution  can  be  ignored  for  kilovolt  electrons  due  to  /, 
being  so  strongly  forward  peaked.  The  quantity  I,1  is  slowly 
decreasing,  as  can  be  demonstrated  by  noting  the  form  of 
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(sin2  By,.  For  the  given  /„  it  can  be  shown  to  be 

(sin2  8),  -  4tjX|lfl  I),!  -  2)  (A20) 

Since  n,E,  is  constant,  the  first  part  of  (AI9),  i.e.,  E/Eo( sin2  8), 
varies  as  /In  q,\  -  2.  This  is  a  decreasing  function  of/,  since  r/, 
is  increasing  from  a  starting  value  of  about  10'*.  It  is  clear 
that  1 1  +  cos2  &),  decreases  from  a  starting  value  of  2.  Its  ac¬ 
tual  variation  for  the  given  p,  is 

( 1  +  cos2  8),  *  |4  +  2 

The  explicit  form  of  IJ  based  on  the  above  considerations  is 
C  **  4t)„(|ln  i)„(  -  2)  (A22) 

Appendix  8:  Derivation  of  the  Average 
Energy  Loss  per  Collision  1£E 

In  this  appendix,  we  given  the  forms  of  EE,  and  its  averaged 
value  A£  It  will  be  shown  that  EE,  is  slowly  varying,  which 
allows  it  to  be  replaced  by  ZlE  in  (5)  without  introducing  sig¬ 
nificant  error.  EE,  is  defined  as  the  energy  loss  per  collision, 
i.e..  loss  over  an  elastic  mean  free  path  at  the  rth  scattering.  In 
the  continuous  energy  loss  approximation  applied  to  this  work 
its  form  is 


n  {PACOS i>  1  (A21) 


E, 

1 


(Bl) 


which  comes  from  ( dE/dx)Es„  where  dE/dx,  the  stopping 
power,  is  given  by  the  Bethe  formula, 


dE  dimZ  <  E 

~  m  — r e  •*>  t 

dx  mv1  l 


(82) 


and  Es„  the  mean  free  path  between  scatterings,  is 


pVi), 


<‘3> 

The  energy  parameter  ]  appears  explicitly  in  our  final  ex¬ 


pression  for  (r)D  (equation  (8))  and  has  the  form 

1  -  9.12(1  +  1.9>rJ/>)  (B4) 

[Segre,  1964).  We  observe  that  EE,  has  only  logarithmetic  en¬ 
ergy  dependence,  since  E,%  -  Es  (see  Appendix  A).  To  obtain 
EE,  we  may  thus  apply  the  continuum  limit,  giving 

M  T'j  ^  /  ~1|  (B5) 
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